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NOMENCLATURE 


A vector is denoted with an over-arrow (e.g., the position vector r ab ), and the column matrix rep- 
resentation of a vector, or any other type of column matrix, is represented with an underbar (e.g., u). 
Points in the system are designated by lowercase letters, so the position vector r ab in the above example 
is read “the position from point a to point b.” Coordinate frames are designated by capital letters (e.g., 
x A is the x unit vector in the A system). Coordinate-frame rotation rates are represented as ul^yj, 
which is read as “the rotation rate of frame F relative to frame I.” Rotation vectors are always resolved 
in the frame of the first subscript. Nondimensional quantities are represented with an overbar and are 
nondimensionalized by the blade mass, rotor radius, and nominal rotor rotational rate.The sine and 
cosine functions are abbreviated to Sq and Cg, where the subscript is the argument. In the following 
list, the abbreviation DOF is used for degrees of freedom. 
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tail-rotor flap inertia about flapping hinge 

blade inertia matrix 

fuselage inertia matrix 
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first-derivative rotor DOF indices (number selected) 
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tail-rotor flap-spring constant 

mount spring constants 

hinge spring constants 

tail-rotor tangent of <5 3 

inflow matrices 

blade mass 

tail-rotor blade mass 

fuselage mass 

rotor-hub mass 

blade-section moment function constants 
number of blade accelerometer measurements 
number of fuselage accelerometer measurements 

number of blades 

number of blades in tail rotor 

number of fixed DOF 

number of harmonics in periodic system linear coefficient matrices 
number of harmonics in trim blade equilibrium solution 
number of rotor blade DOF 
number of empennage surfaces 
number of fixed inputs 

number of controls in blocks A, C, P, S, and G 
number of rotor outputs 
number of states in blocks A, C, P, S, and G 
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SUMMARY 


A mathematical model of a helicopter system with a single main rotor that includes rigid, hinge- 
restrained rotor blades with flap, lag, and torsion degrees of freedom is described. The model allows 
several hinge sequences and two offsets in the hinges. Quasi-steady Greenberg theory is used to 
calculate the blade-section aerodynamic forces, and inflow effects are accounted for by using a three- 
state nonlinear dynamic inflow model. The motion of the rigid fuselage is defined by six degrees of 
freedom, and an optional rotor rpm degree of freedom is available. Empennage surfaces and the tail 
rotor are modeled, and the effect of main-rotor downwash on these elements is included. Model trim, 
linearization, and time-integration operations are described and can be applied to a subset of the model 
in the rotating or nonrotating coordinate frame. A preliminary validation of the model is made by 
comparing its results with those of other analytical and experimental studies. This publication presents 
the results of research completed in November 1989. 

INTRODUCTION 

Future requirements for helicopter flight performance place greater emphasis on high-bandwidth 
flight-control systems than current requirements. Low-order rigid-body models of helicopter systems can 
be used successfully to design feedback control systems, but these models impose bandwidth limitations 
on the controller, making them unsuitable for the design and evaluation of high-bandwidth systems 
(ref. 1). To meet the new requirements, higher order helicopter models must be used in the design of 
the flight-control systems. The higher order effects of primary interest are those of the rotor dynamics 
and the rotor inflow. In addition to capturing the blade-flap and lag dynamics, the effect of torsional 
dynamics on the control design problem needs to be understood. This last item is relevent to rotor 
designs using hingeless and bearingless rotors and to advanced concepts, such as servo— flap control, 
which require torsionally soft rotor blades (ref. 2). 

This report documents a higher order helicopter mathematical model created to fill these needs. The 
helicopter model is described in detail in the next section. The rotor consists of rigid blades with flap, 
lag, torsion, and pitch motions; and with linear hinge springs and dampers. This gives the model the 
rotor dynamics of interest and allows for approximate modeling of hingeless rotor systems. The flap (f), 
lag (1), torsion (t), and pitch (p) hinge motions can have the order Apt, fplt, pflt, lfpt, lpft, or plft (only 
Apt and lfpt are currently available), and two hinge offsets are available to model various articulated 
configurations. The blade-section aerodynamics are linear, and the quasi-steady Greenberg model is used 
(ref. 3). Unsteady inflow effects are included using the three-state nonlinear Pitt/Peters dynamic inflow 
model (ref. 4). An rpm degree of freedom (DOF) of the main rotor is also available, making analysis 
of engine/rotor interaction possible. The rigid fuselage model has six degrees of freedom and includes 
an equivalent-drag-area aerodynamic model, a tail-rotor model, and an empennage-surface model. The 
empennage-surface models have main-rotor downwash effects that are calculated from the flat-wake 
model described in chapter 2 of reference 5. The tail rotor is the quasi-static flapping model described 
in reference 6. For greater flexibility, the spring/damper model, the blade-section aerodynamic model, 
and the fuselage aerodynamic model can be replaced with user-defined models. 

Following the description of the model is the Operations section, which details the operations that 
are performed on the model. There are three operations: trim, linearization, and time-integration. Trim 



places the helicopter system in some user-defined flight condition by satisfying the system equations of 
motion. After trimming, the model may be linearized about the trim condition, or the model may be 
initialized in the trim condition and then time-integrated. Linearization uses nonrotating rotor coordinates 
through a multiblade coordinate transformation (MBCT) (ref. 7). In forward flight, the periodic system 
matrices can be extracted from this linearization process. The time-integration uses rotating blade 
coordinates and allows for arbitrary inputs to the model controls. Additionally, the integration process 
has been cast in a general block form, allowing the introduction of user-specified actuator, sensor, 
feedback, and precompensator dynamics to the system. 

In the Results section, a preliminary validation of the helicopter model is made by comparing 
results with those from other analytical and experimental studies. The rotor model results are compared 
with trimmed rotor-response results, from reference 8 , of an isolated rotor. Rotor/fuselage results are 
compared with the experimental frequency and damping data presented in reference 9. The full helicopter 
model is validated by comparing trim and time-history data extracted from flight data on a UH-60A 
helicopter (ref. 10). Additionally, a comparison with the quasi-static eight-state linearized model from 
the GEN HEL program (refs. 11,12) is made. 

HELICOPTER EQUATIONS OF MOTION 
Fuselage 

The fuselage is modeled as a rigid body whose motion is defined by the six degrees of freedom 
shown in figure 1. The fuselage coordinate frame F is body fixed, with its origin at point m. The 
location of point m is chosen by the user. The body translational motion is defined by the components 
of the inertial position vector rj m , and body rotation is defined by the Euler angles. Other relevant points 
on the body are shown in figure 2. The positions of these points are chosen by the user in the body-fixed 
coordinate frame defined by the station, the buttline, and the waterline. The forces and moments at 
m result from mounting restraints, which can be set to simulate any mounting condition, such as in a 
wind-tunnel test. Point c is the center of mass of the fuselage, where inertial and gravitational forces 
are resolved, and point a is where fuselage aerodynamic forces are applied. The n s empennage surfaces 
produce resultant forces at the Sj points, and the resultant tail-rotor force is at point r. There are . 
points designated as y n , where accelerations are output measurements. Point h is the attachment point 
of the main-rotor hub, which can be set at tilt angles about point t. These tilt angles are defined by 
the transformation, shown in figure 3, which tips the hub coordinate frame H by constant pitch and roll 
angles, 6 S and ® s . 

The equations of motion for the fuselage are formed by summing the forces created by the main 
rotor, fuselage inertia and gravity, fuselage aerodynamics, the fuselage mount, the tail rotor, and the tail 
surfaces. This gives the vector equation 

n s 

Qh + Qig + Qa + Qm + Qr + 5Z Qsj — o ( i ) 

j 

Summing the respective moments about point m gives 

Lfj + Ljq + L i 4 + Lm + Lr 4- Lg =0 (2) 
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where 


Ljj = Jh + ifmt + r t h) x Qh 

(3) 

Ljq = Jj + ( Tmc x Qig) 

(4) 

La = J A 4" (j'ma X Qa ) 

(5) 

Lr = Jr + (?mr x Qr^J 

(6) 

n s 

Ls = 53^ + msj x QSj) 

(7) 
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The final fuselage equations are resolved such that they correspond to the fuselage generalized coordi- 
nates. That is, if a Lagrangian approach had been used, the same fuselage equations would have resulted. 
This resolution of the equations facilitates the extraction of the system mass matrix, which is discussed 
in a subsequent section. The definition of the equations and their location in the fuselage-equation 
vector f are 


feq(l) = ^moments about x F 

(8) 

feq{ 2) = S moments about Ve y 

(9) 

f eq ( 3) = 53 moments about zg z 

(10) 

/eg (4) = 53 forces in x/ direction 

(11) 

/eg (5) = 53 forces in yj direction 

(12) 

/ e g( 6) = 53 forces in zj direction 

(13) 


Using the position vectors and coordinate frames as defined, the velocities and accelerations can 
be determined. At the main rotor hub, the angular velocity and acceleration are given by 


UR/! = & H/F + &F/I 
u F /i = u'f/I 

£ Hj I ~ u'h/I = u 'h/F + Uh/I x ^H/F + & F/I 


(14) 

(15) 

(16) 


Taking advantage of the constant components in the position vectors, the velocities at selected points 
are 


Tit = Tim + Qpji X r m t 
r ih = ? it + &H/I x r t h 
r isj = + &F/I x ^msj 


(17) 

(18) 
(19) 
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r ir ~ r im + ^Ffl x r mr (20) 

r ia ~ r im + &F/I x T ma (21) 

r ic = r im “1“ ^F/I x r me (22) 

r iy n ~ U F/I x r my (23) 

and the acclerations at selected points become 

nt = rim + u'p/j X f mt + Qp/j X (u F/J x f mt ) (24) 

?ih = ru + u' h/i x r th + u H /j x (0 H /j x f th ) (25) 

r ic = r im + ^ Fjl x Fmc + & pjl x ( w F /l x F mc ) (26) 

Fiy n = Fim + u> F/I x Fmy n + ^F/I x i^F/I x 'mt/ n ) (27) 


The forces and moments of the fuselage are discussed separately, beginning with the inertial and 
gravity effects. At the fuselage center of gravity, the inertial force and moment vectors are 

QlG = Mf (-Fic + 9Zi) (28) 

■?I = -Ip ■ uj' pji - £> f /i x (I F ■ up/i) (29) 

where the inertia tensor is defined in the F system. 


The fuselage aerodynamic effects are based in the wind-axis coordinate frame W defined in figure 4. 
The velocity of the air relative to point a is 

= — r ia (30) 

which also defines the direction of the x\y component of the wind-axis frame. The angles a j and 
are defined from the relative velocity as 


ctf = tan 



— 7T < a p < 7T 


(3 f = tan 1 



The fuselage aerodynamic force and moment are 


(31) 


(32) 


Qa — Dx\y + Y y w + Lz w = -pill-alps ’ j (Cp, f x w + Cy f yw +Cp f zw) (33) 

Ja — -Rx\V + Mihv ~ {—CRfXW + ^MjVW - (34) 


4 



which can be converted to the F frame via the transformation Tp\y, which is defined in figure 4. The 
force and moment coefficients are functions of otj and /3y, but the default fuselage-aerodynamics sets 
all force coefficients to zero except Cp } , which is set to unity, making the fuselage reference area Sj 
the equivalent drag area. 

The fuselage mount moments are modeled as linear nonorthogonal springs, and the mount forces 
are treated as given quantities. These forces can be used during trim to satisfy a constraint condition, 
such as setting a certain amount of thrust. The parameters defining the mount effects can be set to zero 
for free-flight analysis. The force and moment expressions for the mount are given by 

L m = (M x - M z Se y ) x F + ( M y Cg x + M z Sg x C 0y ) yp + (- M v Sq x + M z Cq x Cq v ) z f (35) 
where 

M x = -K Fx (6 X - e XQ ) - C F J X (36) 

My = —Kp y (0 y — #yo ) c F y Qy @7) 

M z = —Kp z (9 Z - 9 Z0 ) - C F J Z (38) 

Qm = FxXi + Fyyi + F z zj (39) 

The measurement outputs from the fuselage are acceleration, velocity, and angular rates. The n.Af us 
accelerometer outputs are expressed as 


VU) = (% ~ 9*l) ■ J = l >2,...,n Afua (40) 

fj = fxj xp + f yj yp + fz 3 ZF (4 1 ) 

where the direction of the measurement is selected through the components of fj. The velocity of the 
fuselage center of mass is resolved in the fuselage fixed frame and is available as a measurement, as 
are the fuselage angular rates. The components of the velocity and angular rate vectors are stored in y 
according to 


y{ n Afus + 0 =f ic- *F 

(42) 

y { n A fus + 2 )= f ic ■ VF 

(43) 

y ( n A fu9 + 3 ) = Tic • *F 

(44) 

y{ n A fus +4) =“Fir *F 

(45) 

y ( n A fus + 5 ) = &FH • yp 

(46) 

y ( n A fus + 6 )= $ f/i ■ *F 

(47) 
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Main Rotor 


Fixed to the main rotor are n ^ coordinate frames, one for each blade. Figure 5 shows the rotor 
coordinate frame R for the fcth blade, in which several additional blade coordinate frames are defined. 
The additional frames are shown in figure 6 . The first hinge is offset a distance e from the hub center, 
and point e corresponds to the origin of the link coordinate frame L. The second hinge is offset a 
distance / from point e. The origin of blade coordinate frame B is at point / and is fixed to the blade. 
The significant points on the blade are b, at the blade center of gravity, and x, at a distance x along 
the xg axis. Also on the blade are r>A hld points designated tv n , which locate the blade accelerometer 
measurements. 

The hinge sequences that define the orientation of the blade coordinate frames can be divided into 
two groups . The first group is flap-lag-pitch-torsion (Apt), flap-pitch-lag-torsion (fplt), and pitch-flap- 
lag-torsion (pflt), in all of which the flap motion precedes the lag. The second group places lag before 
flap in lag-flap-pitch-torsion (lfpt), lag-pitch-flap-torsion (lpft), and pitch-lag-flap-torsion (plft) hinge 
sequences. Only the Apt and lfpt sequences are available; they are defined in figures 7 and 8 . An 
intermediate coordinate frame P, which is not shown in figure 6 , is included in these figures. Frame 
P allows the introduction of hinge effects such as precone and 63 . The angle ,3p in figures 7 and 8 is 
the precone angle, and ( p is a cant in the flap hinge that geometrically introduces the £3 effect. By 
definition, the variable £3 is not an Euler angle (ref. 13), but it is related to the cant Euler angle by 
tan 63 = — tan ( p cos 3 p for the Apt system and tan 63 = — tan (c, p + Ck) cos 3p f° r the lfpt system. The 

remaining hinge sequences are defined such that for the flap-before-lag sequences, 3k rotates about the 
coordinates axis —yp\ Cfc. about ~i\ and <£*., about xg. For the lag-before-flap sequences, Q. rotates 
about the axis zp; 3^, about —yi', and 6^, about xg. The sequence of rotations of the blade angles is 
summarized in table 1 . The items in the Apt row and R to P column are read as follows: starting in 
frame R, rotate through an angle of ,3p about the negative yp axis; then rotate through an angle of ( p 
about the zp axis, to the P system. Each element of the table is read similarly. 

The forces on the main rotor are shown in figure 9 with only the A*th blade depicted. A rigid body 
representing the rotor hub mass contributes only inertia to the system. The link between points e and / 
is massless, contributing no inertial forces to the system. Each blade contributes inertial, gravitational, 
and aerodynamic forces and moments, which are related to the loads at hinges / and e by 


Pf = - (Pig + Pa) = -Pt 

(48) 

\lf = - (M/g + M a ) = -Mt 

(49) 

Fe=Ff 

(50) 

Me = Mf + r € f X Fj 

(51) 


The vectors Fj and Mj are the sums of the blade inertial, gravitational, and aerodynamic forces and 
moments at point /. The hub forces and moments of equations (1) and (2) for the main rotor are then 

n b 

Qh = Pig + H {Pt) k (52) 

k= 1 
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This gives the rotor equation 


jfc= 1 


feq( 7) =T q + J H • z R = moment about z R 

(54) 

where T q is the main-rotor torque. For the Apt, fplt, and pflt systems, the flap, lag, 

and torsion equations 

are 

r e q(k, 1) = - M e + Mp + r e f X Fj] • yp 

(55) 

r eq (k,2) = Mf + Mp] • *L 

(56) 

r eq (k, 3) = Mf + Mp] • xr 

(57) 

For the lfpt, Ipft, and plft systems, these equations are 

req{k, 1) = “ Mf + Mp) • VL 

(58) 

r eq (k, 2) = M e + Mp + r e/ X Fp] ■ z P 

(59) 

r eq (k, 3) = Mf + Mp] • xr 

(60) 


These rotor blade equations have been written such that they are equal to the /?£, C fc* and equations 
that would result from using a Lagrangian approach. The total aerodynamic forces and moments at the 
hub center are 


n b 

Sf a = £ (h) t 

k— 1 

n b 

Sm a = £ Wa + { r he + r e f) X F A ] k 

k= 1 

The angular velocities and accelerations are 

ur/i = ur/H +uh/I 

up/I ~ up/R + u R j i 

u'l/I = u'l/p + up/i 

u'B/1 = u'b/L +ul/I 

UR/I = u'r/i = u'r/h + Ur/ I x ur/ H + u H/1 
Up/I - u' P/I = u'p/R + up/ 1 x Up/R + Ur/! 
U'L/I = u'l/I = u'l/P + U'L/I x u L /P + Up/j 
u B/I = u' B/I ~ u'b/L +UB/I x UB/L+ U'L/I 


(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 



For the blade aerodynamic calculations, the velocity of point x results from 


(71) 

(72) 

(73) 


Ae = Ah + w/j// x f he 
Af = r ie +u L/I x r ef 


Ax ~ ?if + v B /i X ff x 

For the mass matrix calculation, the velocity of the blade center of mass is 

Ab = ?if + /i x r/6 (74) 

The accelerations of the blade center of mass, the points w n , and the point x along x B result from 

Ae = f ih + u' RjI x f he + u R/I x (u R/I x f he ) (75) 

r if — r ie + u ' L/I x Af + V L /i X (p R /I x Af) (76) 

fib = Af + u'q/j X ff b + Qq/j X (pBfi X ff b ) (77) 

= Tt/ + u> B/I x ^/u; + U B /I X (p B /j X ff w ) (78) 

Ax = H/ + B/I x #7* + ^S// x ipB/I x ^fx) (79) 

The rotor hub center of mass is at point h, and the rotor hub orientation is referenced to the k = 1 
rotor coordinate system. The inertial and gravitational effect of the rotor hub is 

Pig = m r (~Ah + 9*i) ( 80 ) 

and the inertial moment is 

Pi = ~ [Li • u'r/i + Qrjj x (7 r - Wj?//)] (81) 

Likewise, for each blade, the inertial and gravitational force and moment are 

Pig = Mb (~Ab + 9 * 1 ) (82) 

Mig = - 7 b • u'b/i + £b/i x (7b ■ ^b/i)) + Pfb x Pig (83) 

where the moment is about point /. 


Blade hinge-spring restraining moments are of two types. The first type is nonorthogonal, and for 
the Apt, fplt, and pflt sequences the moments are 

— M e • yp = —Ky0 — Cy0 k 

Mf • = ~K z C,k — C z <tk 

Mf ■ x B = -K x <p k - C x <p k 
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(84) 

(85) 

( 86 ) 



For the Ifpt, lpft, and plft sequences, the nonorthogonal moments are 


(87) 

( 88 ) 

(89) 


— Mf • VL — ~ky,dk ~ Cyfik 
M e ■ zp = -K-C k - CzCk 
Mf • ip = -K x <t>k ~ C x <f>k 


The second type of spring restraint is orthogonal and is used to approximate hingeless rotor-blade 
systems. These orthogonal moments are valid only for (1) no second offset (/ = 0), (2) no cant angle 
(£ p = 0), and (3) with the pitch inboard or outboard of the flap and lag degrees of freedom. For the 
Apt and pflt sequences, the moments are 


Mf = Me = 


-K* + 4>k^Ck C 0k) ~ Cx {~^ S Pk + tkQkCpk)] X P 

+ -Ay (-/3fc + <f>kS - Cy (~fik + VP 

+ [~Ar (c kCflk + - Cz {CkCfik + hC Ci S fik )] Zp 


(90) 


For the lfpt and plft sequences, the moments are 


M f = M e 


-K x + 4> k C0 k C Ck ) - C x (0 k S Ck + ik C0 k C Ck )] £ P 

+ [-A y (-PkC<; k + faCphStk) ~ Cy {~Pk C ( c* + 


+ 


-Kz (Ci fc + 4>kS0 k ) - Cz (Cj t + n%)] ip 


yp 


(91) 


The total blade aerodynamic effects are found by integrating the sectional aerodynamic forces and 
moments along the span according to 


where 


- f x 2 - 

Fa = Sa dx 

Jx\ 

M\ = J x " [(* ip) x Ia + 


dx 


Ja ~ fyc^C fzc x C 
in A = m x c ic 


(92) 

(93) 


(94) 

(95) 


The integrations are carried out using Gaussian quadrature numerical integration. Blade-section forces 
are resolved in a chord coordinate frame C shown in figure 10, which is rotated from the blade coordinate 
frame B by a geometric twist. The linear twist function is defined by 

@t( x ) = TwxX + T wc (96) 
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The components of the sectional forces and moments are expressed in terms of the angular velocity and 
acceleration of the chord section and of the components of the relative velocity and acceleration of the 
air mass at the origin of the C frame. The latter two are 

Vc = V G + A,- - f ix (97) 

Ac = Vc~ n x (98) 

where the vectorA; is the induced air velocity whose derivative has been ignored in the air-mass ac- 
celeration. The induced air velocity results from the dynamic inflow model, which is described in a 
subsequent section. 

The sectional aerodynamic model is based on quasi-steady Greenberg theory, which is a Theodorsen 
theory modified to account for lead-lag motions (refs. 3,14). The lift and moment about the torsion axis 
are given by 


L — PA a b ^ Q 4 " P Aero 

AI = PA a b f J-ab "l” 2 ) ^ Q {x ab P ^ 8 2 ? -^ ero 


where 


Q = (h + Va) - - 

P = (A + Vaj - x ab a 


x ab 2) a 


( 100 ) 


( 101 ) 


( 102 ) 


The variable i Aero 1S a switch to remove the apparent mass effects, which can be ignored when the 
order of magnitude of the mass of the airfoil section is much greater than the order of magnitude of the 
mass of the air it occupies (ref. 14, chap. 5). Drag and additional moments due to camber are added by 


D — dQ + d\ a e ff + dj 
Alcam = m 0 + mi a e ff 


(103) 

(104) 


Figure 1 1 shows the application of the sectional aerodynamic model in forward and reverse flow con- 
ditions. The variables in the Greenberg theory are interpreted as in reference 15, where the quantity 
h + Va is interpreted as the normal air velocity at the torsion axis. The derivative of this quan- 
tity is then the derivative of the z c component of the air velocity V' c , which can be extracted from 
A c = V'c + ^cjl x VC- 

Blade measurements are stored in w in the following manner. For the A:th blade, the n^. 
acclerometer outputs are given by 


[k + n b (n- 1)] = ( f iWn -gzj) 


n = 1,2, ... , n i 


K - b Xn x B + by n y B + b Zjl z B 
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where the direction of the measurement is selected through the vector b n . Blade moments at the hinges 
are also in w . They are given by the following expressions when flap is before lag: 


w 


w 


w [k + n b n A^i ( i 
k + n b (n AbU + l) 


k + n b ( n A bld +2 )] 
and by the following expressions when lag is before 


w 


k + rib n A bld 


w 

w 


k + n b (n Am + l) 
k + n b { n Abld + 2 ) 


X 

+ 

(107) 

M t • Z L 

(108) 

Mp ■ x q 

(109) 

flap: 

M t • yi 

(110) 

[Mj + r e f x Fp^j • zp 

(111) 

M t ■ x B 

(112) 


Inflow 

The rotor inflow is modeled using the three-state nonlinear Pitt/Peters model of reference 4, which 
is applied in the actuator-disk coordinate frame A shown in figure 12. This frame is tipped to the left 
by the sine of the flapping angle and tipped forward by the cosine of the flapping angle. The velocity 
and acceleration of the air relative to the hub center are 


V h = V C -'r ih (113) 

A-h = V G ~ ^ih ( 114 > 

The advance ratio, the free-stream inflow, and the total average inflow can be found by using the velocity 
vectors; they are 


M • *a ) 2 + W h v *) 2 
^ = m 

x YhjJa 

fs m 

^ ^ fs + ^0 

where the total average inflow A is the sum of the free-stream inflow plus the constant induced-inflow 
component, which is one state of the dynamic inflow model. The inflow vector for the blade aerodynamic 
calculation, in terms of the inflow states, is 

A,- = Xiz A = (a 0 4- Aj s — Z A (118) 


(115) 

(116) 
(117) 
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equations 


* > 



v 0 

* 1 

r -° l 

{ C T 1 

^ 1 

r vu 

' ~ L I \ c L 

v l c J 

\ •'lc J 

{ C M J 


(119) 

th V”' l0W in ' he $ame Way that <*«"' the 

velocity in the d ts k p.aL The co Z ^ “ 


= Ao 

Uls = X isC/3 h + Aj c Sp h 

v l c ~ X lcCp k + A \ s Sp h 

Ao 
fi 


* 

^0 


; ‘ s = lf% + yr s h + 


£2 


- f + AjcC^) 


fi v 


1/10 C 0h Q s P h ( A ic‘S/? A + A^C^J 


where the side-slip angle and rate are 

0A = tan - 1 


= 


~ V h • y.4 

“t'A ' XA t 

' Xa4 ) (j* ' $a) - (4 • X A ) (fy, ■ 

(^A ' x a) + (^A • Va) 2 


( 120 ) 

( 121 ) 

( 122 ) 

(123) 

(124) 

(125) 


(126) 


(127) 


The thrust, roll, and pitch coefficients in equation (1 19) are eiven in t * rmc t u 

forces and moments, by ) g , m terms of the main-rotor aerodynamic 


i 

ii 


Kp A R 4 fi 2 


r L _ &U A ’ Va) % + (s Ma ■ x A ) 

1 w 

t rp A R 5 Q 2 

M off 

r u _ ($Ma ' Va) C 0h - (5 Ma ■ x A ) 

S Ph , , 

t, P a R 5ft2 

M off 


(128) 

(129) 

(130) 
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The matrices for the Pitt/Peters model are 


L I = 



0 


0 


-4 

j/( 1+5 0j ) 


15t 

64j/7’ 



2* 


0 


15tt /l , 

o 


-4S Q , 

"(l+So,) 


Mj = 


m.Q(mass ) 0 0 

0 ^ 0 
0 0 ^1. 


where 

__ /i 2 + A (A + Ap) 

v/TTaJ 
= >/?^ 
a s — tan 



(131) 


(132) 


(133) 

(134) 

(135) 


The element mg takes on two possible values. One value corresponds to a zero-inflow boundary 
condition at the hub center, mo(0) = jjf ; the other value, mo(l) = does not enforce this condition. 


Tail Surfaces 

Up to five tail surfaces are available (the subscript representing each surface will be omitted in 
the following discussion); they are oriented as in figure 13. Each surface can be tipped away from the 
horizontal by an angle a x , and the pitch incidence can be adjusted through a y to define the zero point 
of the input ixs- The velocity of the air relative to the surface is 

V s = Vq + drs ~ ?is ( 136 ) 

where djg is the downwash from the main rotor. Each tail surface has its own downwash table that is 
a function of the angles at which the air strikes the hub center, the magnitude_of the air velocity, and 
the thrust coefficient (i.e., d TS = CtSdW ts ( q v, Pv,V v ) z s ). The moment J s at the surface is zero, 
but the force at the tail surface is 

Qs PA II V s || 2 £ts{[Cz, (<*;) - Cd ( a - r ) Ca s \ is 

+ [ — C i (o z ) Ca s ~ C D ( a z) 9/oss *-*/? (137) 
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where, for forward flow, a z and S R are defined by 


and, for reverse flow, they are 


where 


Otz 


a s 

Sr 


1 

defined 

by 


<*z 

= 

a — 7r 

Sr 

= 

0.8 

a z 

= 

a + tt 

Sr 

— 

0.8 


a s 

= tan 


-1 


> for — ^ ot§ ^ 

for 5 < tt 

for —7 r < a s < — ^ 



( 138 ) 


(139) 


The lift and drag coefficients are functions of a z and are defined by specifying the break points shown in 
figures 14 and 15. The dynamic pressure loss due to fuselage blockage is approximated by the functions 


Qloss — 1 


^/oss \ 
Vs ) 


2 


Vs 


= C TS ex p 



f Q / Q 0 TS \ _j_ 

V a TS ) 


{ Pf_ ^0 TS \ 

\ hs ) 


(140) 

(141) 


When the angles aj and j3f are at ao TS and bQ TS , respectively, the loss is a maximum. 


Tail Rotor 

The tail-rotor orientation is shown in figure 16; the Euler angles p x and p y have the same sequence 
as the tail-surface orientation. The relative velocity of the air is given by 


v r = Vq + d TR -f ir (142) 

and the downwash dj R , like the tail surface downwash, is interpolated from a table. From the figure, 
the advance ratio and total inflow are 


TTR = 


(Vr ■ xr ) 1 + (v r ■ y T y 


Rtr^tr 


a TR = VTR tan a r + 


Cj. 


TR 


2'J^TR + 


“ X f»TR + 


Cj, 


TR 


TR 


2yTTR + ^ tr 


The wind orientation angles with respect to the hub are defined as 


jS r = tan 1 


\-Vr- &T) 


(143) 

(144) 


(145) 
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(146) 


a r = tan 


-1 


[y(f r .x r ) 2 +(F r .y t)‘ 

The force and moment for the hub from the tail rotor are 


V r • z T 


Qr = ( ~ H TR^p r ~ Y TR^p r ) XT + {~ H TR^p r + YTR C P r ) VT ~ T '. TRfblk*T ( 147 ) 

Jr = (j-TR^pr ~ M TR S p r ) XT + ( L TR S p r + M TR C p r ) VT + QtR~ z T ( 148 ) 

All components of the force and moment expressions result from the quasi-static flapping model of 
reference 6, and are functions of the total inflow and the blade harmonic components. The inflow 
comes from equation (144) and the blade harmonics result from 

[K}a-f = 0 (149) 

where a is a vector containing the collective, the first sine, and the first cosine components of the blade 
flapping motion. The precise definitions of the matrices and the rotor forces and moments are given in 
appendix C of reference 6. In the present report the total inflow is defined such that positive inflow 
is down through the rotor, whereas in reference 6 total inflow is defined such that positive inflow is 
up through the rotor. Both equations (144) and (149) are algebraic nonlinear equations and are solved 
using the Newton method. Once acceptable convergence is reached, the values of a and Ajtj can be 
used in the force and moment equations. The loss factor /&/£ in equation (147) is included in order to 
model the interference effect between the rotor and the surface to which it is mounted. Loss factor /&/*. 
is defined as follows: 


fblk — (TRloss 0 
fblk = 1 


I 1 {trL) +1 


PTR < TRbrk 
PTR > T R brk 


(150) 

(151) 


The variable TRi oss is the percentage of thrust available at zero advance ratio, and TR^j. is the advance 
ratio at which / ^ becomes unity and remains unity at advance ratios above 


Main-Rotor Downwash 

The downwash of the main rotor on the tail surfaces and tail rotor is calculated using the flat-wake 
model described in detail in chapter 2 of reference 5. The wake model assumes a rotor in edgewise 
flow, as shown in figure 17; the coordinate frame is that of reference 5. The downwash velocity is 
determined using the Biot-Savart law, assuming cycloidal vortex filaments shedding from the disk plane. 
This flat-wake assumption is considered valid in forward flight when VjRVt > \.63y/Uj The x, y, 
and 2 components of the velocity are calculated from 

^(xi,yi,Ii) = ^=^ 2 A£ )> ^dp j = x,y,z (152) 
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where Avj results from 


AO* 

Av y 


Av z 


(x h yi,zi) 

A v yx (*l,yi,«l) 


- ^ ^ + (oo»yi> 5 l)] + [iV(5i,yi,zi) + iV(oo,yi,2i) 

[A' (xi,yi,fi) + K (oo,y 1( 5i)] + [£ (ii,yi,5i) + l(oo,yi,£i) 


(153) 

(154) 

(155) 


The functions in equations (153) through (155) are pretabulated and depend on the nondimensional 
distance from the hub center p and the components of the point of interest, through 


*1 


x \ • ,-7 Vl 

— ; yi = — 
p P 




fi 


(156) 


A typical example of each of these functions is shown in figure 18 for y\ = 0.6. The functions need to 
be tabulated only at the positive values of x\, y\, and z\ because of the symmetric and antisymmetric 
properties of the functions. A summary of the odd and even behavior of the functions is shown in 
table 2 along with the equation numbers and figure numbers that define the functions in reference 5. 
The downwash components are expressed in terms of the function T, which is the average circulation 
over a rotor revolution at a radial distance from the center of the rotor. As in appendix C of reference 1 6, 
this function is approximated as 


r ( ?) = r 0 ( l - 7J)p 2 



20ttC t 

5 ( r 4 - r?) - 4 (r| - rf) 


(157) 

(158) 


where the constant Tq is a function of the thrust coefficient. This relationship between the circulation 
and the thrust coefficient comes from evaluating 


T = 


r r 2 ' 1 & 1 
Jr\ 2 tt JO 


2 1 dr 

dp 


dtp 


dp 


where _ 

is: 

which results from the Joukowsky lifting law. 


(159) 


(160) 


The wake model can be used to generate a downwash table for each tail surface and tail rotor. The 
wake is oriented along the relative air velocity vector in the xy — yy plane, shown in figure 19. The 
S-frame is the tail-surface frame, but the tail-rotor frame T can be substituted for it. The velocity at the 
hub is 


V v = Vq + \qza 


~ Uh 


(161) 
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and the orientation angles are 


a v = tan 



0 V — tan 


( -Vy • m ) 

\-^u • % / 


( 162 ) 


(163) 


For the given tail surface (tail rotor) fixed in the S (T) frame, downwash values fow TS Ud\V T r) at 
different combinations of a v ,/3 v , and ||K,|| are calculated and stored in a table. These values are the 
downwash velocities in the z§ (zp) direction divided by Cp . This division of the downwash velocities 
by Cp is based on equations (152) through (158), from which the proportionality of C 'p to the downwash 
velocity can be deduced. For the tail rotor, the downwash vector is 

dpp = CTfDW TR { a v,0v, tlKill) *T 

and for each tail surface the downwash vector is 

dp$ = CtIdWts ( a v,fiv, ||Vu||) % (^5) 

The fow functions are linearly interpolated from values tabulated from a given set of a v , (3 V , and V v . 

The points that make up the interpolation tables are selected in the following manner. The ||V u ||’s 
are chosen directly by selecting the number n VDW and the particular values. The q^’s are chosen by 
selecting a set of points a pp and defining = a. pp + ct v q; similarly, the (5 V s are chosen by selecting 
a set of points j3pp and defining f3 v — 3pp + 3u0- The quantities q^o and (3 v q are the orientation angles 
that center the wake on the given surface (i.e., when V v is parallel to r ). The points app are chosen 
by setting the boundary angles oliq bnd < a iq < — &DW < 0 < &DW < a UP < Pbnd' an< ^ 
points j3pp are chosen by setting the boundary angles 0lo B nd < @ LO < ® < @UP < @UPbnd' These 
boundary angles are defined for the condition in which the flat wake is centered on the given surface. For 
the a boundaries, the number of divisions in the end regions, [■ aio BND ^ a LO ] and i a UP’ a UP B ^ D ^' 
is selected through n aend , and the number of divisions in the middle regions, [a L0 , ~&Dw] and 
[A DWi a UP\' * s selected through n amid • The value A p)\y is a small value and precludes evaluation 
of the downwash on the wake layer, where the vortex singularities are present. For the 3 boundaries, 
the number of divisions in the end regions, [{3lo bnd , 9 lo\ and [&UP-> &UP B hd\' IS se l ected through 
no , and the number of divisions in the middle regions, [/?£o,0] and [O./W], is selected through 

° end 
n Pmid- 


Mass Matrix 


The equations of motion are formed numerically, making the model described in this report implicit. 
Because of this, the system mass matrix is not readily available as it would be if the equations were 
derived explicitly. Allowing for different hinge sequences precludes the explicit calculation of the mass 
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matrix, since each of the many hinge sequences would require a separate calculation. To avoid this, the 
mass matrix is implicitly calculated with the condition that no acceleration-dependent applied forces, 
such as the apparent mass terms in blade-section aerodynamics, are present in the system. Such forces 
contribute to the mass matrix, making its extraction more difficult. It is well known from analytical 
dynamics that the form of the kinetic energy is 


^kin ~ ^2 "T T\ + Tq — q fqq(q, t) q + C(q , t)^q + r(q, t) 


(166) 


where 


the model presented. 


Mqq is the mass matrix of the system and q represents the system generalized coordinates. In 


T 

Q = Pi, P2) ••• >Pn b , Cl >C2» ••• i Cn b ,01,02> ••• , Rm x 3 R m y 3 Rm z , 

' v 

rj jT 

and the quadratic contribution to the kinetic energy is given by 

n b 

^2 ~ ^2 f use i age + ^2 rotor + X/ (' ^blade)k 

k=\ 

For the fuselage, the velocity and angular rate are given by 

W F 


= W F/I fi = 


U F/I '■ 


:0 


0 

R 


.(F) a , .(F) „ A L(F) 

T-ic ~ T FlE + r mc Wp/j© = rjnc Uf/i ■ Tpj 


0 

i? 


(167) 


(168) 


(169) 

(170) 


R f 


which is a re-expression of equation (22) using matrix notation instead of vector notation. The elements 
of the column matrix R are the time -derivatives of the three translational degrees of freedom describing 
the motion of point m on the fuselage; the elements of 0 are the time-derivatives of the three Euler 
angles of the fuselage. For the coordinate system F defining the fuselage motion. 


U F/I = 


1 0 

0 C 9x S 9x c 0l 

0 -s 9x c e c 9 


(171) 


The velocity and angular rate give the T 2 contribution to the kinetic energy from the fuselage 

T 


T 2 fuselage = \ { | f [M F R T f R F + W T p I F W F 


1 } 


(172) 
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Likewise, for the rotor hub mass, the matrix representation of the velocity and angular rate is 


a = s±4y / + y^H/F — R/H ~ TrF^F/I® + & R/F 4> ~ ^RF^'F/I '• 0 : “-’fl/fj j 4 } (173) 


r -,P = T F/I& + ^mh^F/lQ - f mh^F/I : T F/I : 0 \ 4 


(174) 


where 


—R/F ~ 0 4'— Ur/ pi' 

1 


(175) 


Equations (173) and (174) arise from equations (14), (18), and (63). As with the fuselage, the T 2 
contribution to the kinetic energy from the rotor hub is 


T 2rotor = ^4 [MrR T r Rr + W T R I R W R \ 4 


(176) 


The same approach is used to find the kinetic energy for the fcth blade. For the fcth blade the angular 
rate is 

= P R /\ + ^l/r + ^-B/L ~ TbF&F/I® + T BR u R /f4' + T BL u L / R a L / R + u B j L a B j L 


= \Tspupjj : 0 : T BR u R /p : TbL&l/R, : &B/L 


i IL/R 
^B/L 


and from equation ( 74 ) the velocity of the center of mass becomes 

(tf) (fl) (fl) 

Ue ~ Uh + L he y-R/I 

AL)_AL) AD 

^ if ~~ —it ' —ef —L/ 1 


+ ^fi^-B/l = MyUF/i : T B j : Mju^p : M\U L j R : rjpj ) u B / L < . ij> 

' ^ L/R 

RBk l ^B/L 
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where 


Mi=T BL rl?+r ( f B b ) T BL 
M 2 = M\T lr + T BR f[^ 
M's = M 2 T RF + T BF f ^ 


Since no particular hinge sequence has been stipulated, the (Lij R and a B j F are generic, representing 
the generalized coordinates in the transformation indicated in the subscripts. For the Apt, fplt, and pflt 
hinge sequences, the terms are defined as 


^-L/R ~ 


—B/L = 


5 C p 


-a 


c p 


0 


0 


S 9 k +4> k 

Cd k +<p k 


h = “L/R^L/R 



“B/L^B/R 


and for the lfpt, lpft, and plft hinge sequences, the terms are defined as 


^-L/R = 


' 0 
0 
1 


C k ~ ^L/R^L/R 


^-B/L 


Cd k +<p k C<; p 0 
^9 k +<t> k C^ p 0 


Pk \ _ - 

k j~ UJL f R - B / R 


(179) 


(180) 


(181) 


(182) 


Using the velocity and angular rate, the T 2 contribution to the kinetic energy from the kth blade becomes 



Q. 1 

T 

' 0 ' 

R 


R 

0 

• [M B R T B k R Bk + w Bt i B w Bt ] , 

■0 > 

&L/R 


—LjR 

k ^b/l t 


. 9 b/l t 


(183) 


For the inflow model, the mass matrix associated with the first-order states x T = [Aq, Ai s , Ai c ] can be 
easily extracted using equation (119) and equations (123) through (125). The mass matrix for the inflow 
model is 


Mdi 


1 

n 


LiM! 


10 0‘ 

0 C 3 h % 

0 -% c Pk J 


(184) 
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If there are no acceleration-dependent applied forces, the total mass matrix of the system can then be 
related to the equations of motion by 


Mqq 0 
0 Mqi 



(185) 


The matrix M qq is formed using the T 2 terms in equation (168), by adding the contributions of the 
corresponding fuselage, rotor hub, and blade elements. The equations included in are the rotor-blade 
equations, equations (55) through (57) or equations (58) through (60); the equations included in are 
the fuselage equations, equations (8) through (13) and equation (54). Finally, the equations included in 
x_ eq are the inflow equations, equations (119). The reader will note that these equations have been cast 
such that they are identical to equations that would result from a Lagrangian formulation. Thus, the 
mass matrix calculation is compatible with the equations of motion and allows the relationships 


-“f = 1 

' Leq 

f 

±-eq J 

' Mdi = M 

\£-eq 




Leq 

[o, 0 , q , x, q ; 

1 

g = < 

f 

. 1 

JO, 0, £, x, £ 

(o, 0, q , x, q t 

! J 


Multiblade Coordinate Transformation 


(186) 

(187) 

(188) 


For rotor systems with three or more blades, a multiblade coordinate transformation (MBCT) is 
available to transform the rotor rotating coordinates to nonrotating coordinates (ref. 7). For a set of 
rotating variables, each of which represents the same physical quantity for each blade, 

3& = [ x U x 2i x 3i"-] ( 189) 

The nonrotating representation for an even number of blades is 

— [-To? ^ Is ’ lei • • •] (190) 

The generic variable xq is the collective, while xj is the differential, which is not present if the number 
of blades is odd. The x ns terms are the nth sine variables, while the x nc terms are the nth cosine 
variables. An example is the flap degree of freedom, for which xr is a vector whose components are 
flapping angles. The transformation between variables is linear and time-varying, so the transformation 
from rotating to nonrotating systems is given by 


(191) 

(192) 

(193) 


ln = t nrlr 

£jV = TfifRXR + TnrX R 

x N = T nr x r + 2 T nr x r + T nr x r 
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and the transformation from nonrotating to rotating systems is 


£R = T RN£-N 

(194) 

i-R = TrnXn + Trn£.N 

095) 

£r - Trn£n + 2 T NR X N + T NR X N 

096) 


The form of the transformation matrix for models with an even number of blades is 
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The ipk terms are the azimuth positions of the blades: 


1>k = ^1 + (k~ 1)— 

n b 


k=l,2,...,n b 


(199) 


These positions all have the same time derivatives. For models with odd numbers of blades, the second 
column of equation (197) and the second row of equation (198) are not present. 


OPERATIONS 


Basic Setup 

A driving routine performs four basic operations on the helicopter mathematical model: initial- 
ization, trim, linearization, and integration (fig. 20). Initialization sets the basic data that define the 
helicopter configuration and conveys to the driving routine all the data necessary to operate on the 
model. The trim operation sets the helicopter configuration into some prescribed steady condition. A 
harmonic balance technique is used, which simultaneously balances the fuselage forces and moments 
and finds the rotor-blade equilibrium solution. Once the model is trimmed, two possibilities are avail- 
able. One is linearization, in which a two-point difference formula is used to linearize the model about 
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the trim condition. The other possibility is time integration, which gives the system time response, 
initialized in trim condition, for given control inputs. 

Each operation uses a subset of the model, whose general structure is shown in figure 21 and 
represented in equation form as follows: 


E.(X) = 


Y = 


X 1 = 


T-eq 

f 

± -eq 

3ieq 

W 

y 

x 


M (r, 


r 

f 

x 


+ G 6,u) =0 


r T ,f,x T ,r T J T ,x T ,r T ,f T ,e T ,u T 


(200) 


(201) 

( 202 ) 


Variables and controls are represented in a single vector X, which acts as the input to the model, whereas 
the equations and measured quantities are represented in a single vector Fj q — [F?Yj\, which acts 
as the output. Both the inputs and the outputs can be categorized into three types: rotor-blade, fixed, 
and augmented. The rotor-blade inputs/outputs are items associated with the rotating coordinate system, 
such as the flap degrees of freedom or the blade pitch; the fixed inputs/outputs are items associated with 
the nonrotating system, such as fuselage degrees of freedom or fuselage acceleration measurements. The 
augmented inputs/outputs are associated with any additional first-order equations necessary to complete 
the model. The structure in figure 21 allows the use of rotating or nonrotating inputs and outputs. 
The variables are transformed by the multiblade coordinate transformation; the choice of the coordinate 
system depends on the particular operation. This basic setup is general and the operations can work on 
any model set up in this manner. 

For the model described in the Helicopter Equations of Motion section, the rotor-blade degrees of 
freedom in the rotating system are 

r T = ( 203 ) 

When using nonrotating coordinates, the vector r is defined according to equation (190) for each degree 
of freedom. For a model with an even number of blades the rotor degrees of freedom in the nonrotating 
system are 

r T = [fio, 3 d , fiis, fi\c, . . . , Co, C d > Cls, Cl c, • • • , 00. <t>\ s> <J>1 c> ■ • •] ( 204 ) 


For a model with an odd number of blades, the differential coordinates are not present. Using nonrotating 
coordinates for the input variables requires that the flag ISNROT = 1 , which transforms the nonrotating 
coordinates to the rotating system via equations (194) through (196). 


The rotor-blade controls in the rotating system are 

0 T = [0i, 0 2 , •••,*1.02, 02,..-' 


(205) 
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and as with the rotor-blade variables, the function of the flag IINROT is to transform 0 to the rotating 
system if these inputs are given in the nonrotating system. For a model with an even number of blades, 
the rotor-blade controls in the nonrotating system are 


£ — 0Q> 01s>#lc> • • ■ , 0O>0d>01s>01o • • • j0O>0<F01s)01c> • • • 


(206) 


For a model with an odd number of blades, the differential coordinates are not present. The fixed 
variables are the fuselage degrees of freedom and the rotor azimuth angle 


/ — ? Qy i Rm x i ^rn y i Rm z ? 


and the augmented variables are the three dynamic inflow coefficients 
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The fixed controls u are 


— [^Orif v g x i v 9y' v 9z 1 ^TS\ > *TS2’ • • • > ^TS ns ] (209) 

The center block of figure 21 is the model, described in the previous section, that calculates the helicopter 
equations of motion. The equations in r eq are the rotor-blade equations, equations (55) through (57) 
or equations (58) through (60), and the equations in / are the fuselage equations, equations (8) 

— cq 

through (13) and equation (54). The equations in x^ q are the inflow equations, equations (119). The 
vector y is defined by equations (40) through (47), and w is defined by equation (105) and equations (107) 
through (112). The vector w is in the rotating frame but can be transformed to the nonrotating frame 
by using the flag IONROT. 


Initialization 

Initialization performs two basic functions, one of which is to set up any data not given in X 
that define the model. These data are the physical parameters defining the model that was previously 
described; they are summarized in table 3. The main-rotor downwash tables are either generated and 
stored in files or read in from previously generated files. The second function of initialization is to 
transfer basic system data to the driver so that subsequent operations can be carried out on the model. 
These data, summarized in table 4, are primarily size data of the full model. Included in the table are 
the current values of these data for the model described in the Helicopter Equations of Motion section 
of this report. 


Trim 

The trim operation places the helicopter in the desired flight configuration by using a harmonic 
balance technique that is described in reference 8, chapter 3. The method casts the helicopter equations in 
an algebraic form by defining the rotor-blade variables as a Fourier series. The coefficients of this series 
expansion, along with the fixed and augmented state variables and inputs, are adjusted to force the rotor- 
blade, fuselage, and augmented state equations to zero. The rotor-blade equations are also expressed as 
a Fourier series, so that forcing the rotor blade equations to zero means forcing the coefficients of this 
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series expansion to zero. It is assumed that the helicopter main rotor has a constant rotation rate ip = 
and that the rotor-blade controls are expressed in the nonrotating frame (i.e., IINROT = 1). All blades 
are assumed identical, and information from only one blade, expressed in terms of rotating coordinates 
(i.e., ISNROT = 0), is needed in the trim. An EMSL nonlinear algebraic equation solver, ZSPOW, is 
used to find the solution to the following equation (ref. 17): 


F(V) = 0 


( 210 ) 


The vector of inputs, V_, is defined as 
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which is composed of fixed, augmented, and rotor-blade inputs. (Matrices are stacked by columns 
when represented in column vectors.) The blade variables are represented as coefficients of the Fourier 
expansion. 


Uh TRM 

r[k + n b (j - \ )}=r co {\,j)+ £ r c0 (2n,j)cos(nip k )+r co (2n+ \,j)sin(mp k ) (212) 

71=1 

The variable k is the blade number, and j is the blade degree of freedom (e.g., j = 2 is the lead-lag; 
see eq. (203)). For the model of this report, the coefficients are represented in matrix form as 
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The elements of V are adjusted by the equation solver after the user chooses determining indices. For 
example, the choice of n lf = 2 with if( 1) = 2 and if (2) = 3 directs the solver to adjust angles d y and 
6 Z to find the trim. Given the information in V_, the input vector X can be formed and the model can 
be evaluated. The resulting equations can then be used to form F, which is defined as 

r fo K M 1 )] n = 1)2, , (2 n hTRM + 1 ) 

r fo [n,jr(n jr )\ n = 1,2, ... , (2 n hTRM + 1) 

% g** Leq [heq ( n )] dt n= 1 , 2 ,..., n feq 
[*l(n)] ti>i (n)dt n = \,2, . . . ,n i{ 

n 

r fo(lJ) = 2 ^J 0 ^eq t l +n b(j ~ l)]dt 

ft r2*/n 

r /o(2> j) = — y o Leq [1 + n b(j ~ 1)] COS (Qt)dt 

Q rlir/Q 

r /o( 3 >j) = -y 0 Legll +n b {j - l)]sin(Qt)dt (215) 

ft f2ir/il 

r f ~ Jq i-eq [1 n b(J ~ 0] COS {n^ ^^Vtt)dt 
f 1 rln/Q 

r f o(^ n hxRM •?') ~ Jq Leq [1 + n bU ~ 0 ] sinfohj- RM^)dt 

These equation are the Fourier coefficients of the first blade equation, where 

Uh TRM 

r eq [& + nfj(j — 1 )] = Tf 0 (\,j) + ^2 r f 0 (2n, j) cos(nij>k) + r f 0 (2n + l.j) sin(n^) (216) 

n=l 

The choice of rih TRM determines the number of rotor coefficients that are adjusted to force the same 
number of coefficients of the blade equations to zero. The fixed and augmented equations are averaged 
over one rotor revolution in equation (214) to zero only the constant harmonics of these equations. A 
summary of the inputs that control the trim procedure is shown in table 5. Indices and their number (in 
parentheses) are selected to choose the adjustable variables and the equations that are to be satisfied. If 
an adjustable variable is set to a value, the value acts as an initial condition to the solver. If the variable 
is not adjusted by the solver, it remains at the set value throughout the trim. An example of the use of 
this feature would be to trim the model at a constant climb rate. 

Integration 

The time response of the system shown in figure 22 to given inputs u p can be obtained by numerical 
integration. Since it is difficult to generalize helicopter control/actuator models, these systems have 




26 



been broadly defined as first-order equations that can be defined by the user. The P block in figure 22 
represents precompensation, and the A block is any actuator/mixer dynamics; feedback dynamics are 
represented in the C block and sensor dynamics are represented in the S block, which has the restriction 
that the output of S cannot explicitly depend on u s . The helicopter model is represented in block G with 
nonrotating rotor inputs (IINROT = 1) and rotating-blade degrees of freedom and outputs (ISNROT = 0, 
IONROT = 0). By default, P and A are identity feedthrough and C and S are zero, giving the response 
of the open-loop helicopter model G. Any subset of the model represented by equation (200) can be 
integrated, where the state of G is defined by 
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As in the the trim procedure, the selection of the variables and equations is done by choosing indices. 
If the rotor degree of freedom is not selected, a constant rotor rate is used, by default, in the integration 
of the model (i.e., ip = Q). The inputs are defined as 
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The feedback outputs are 
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The equations integrated are 
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From equation (200) the form of the subsystem that is integrated is 


- Mq (xq,uq) d,Q + G_g (xg,Ug) = EjG (dG^G^G) 
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where 

f £{n + n 6 [7r(m) - 1]} n = 1,2, m = l,2,...,n jV 


dG = j/zy(m) m= l,2,...,n^ i (223) 

, m= 1,2,...,^, 


Two methods are available for extracting the mass matrix Mq. One method uses equation (222) 
directly, in which it is clear that the jth column of Mq is F_q, with d G equal to zero except for the 
jth element set to unity, minus F_q{0,xq,uq). If Mq is of dimension n\j G — n b rij^ + + n^, 

one evaluation of the model is necessary to get Fq(0, xq,uq), and n \/ G additional evaluations of the 
model are necessary to get Mq. This approach always works; however, it requires many evaluations 
of the model and thus is inefficient. The other method is to use the corresponding elements from the 
mass matrix, which is calculated in and passed from the model, to form the submatrix Mq. Once Mq 
is available, the equation 

Mq 1 (xq, UQ ) Gq (xq, uq) 

i[n + n b jf(m ) - 1] n = 1, 2, . . . , n b ; m = 1, 2, . . . , nj. 

f\if(m) m= l,2,..,,n^ 

can be found; this is the expected form of block G in figure 22. 

If the helicopter model is trimmed, the inputs and outputs of G are available at the initial time. The 
initial state x q q , initial input uq q , and initial output i/q q can be used to extract the initial conditions 
of the blocks A, C, P, and S (fig. 22) as follows. It is assumed that = xq = xp = x§ = 0 and 
it is required that n UG = n Up , which gives an equal number of equations and unknowns. With these 
conditions, a nonlinear equation f{v) = 0 can be set up, where / and v are defined as in figure 23. 
The variables on the left-hand side of the figure are adjusted until the variables on the right-hand side 
are forced to zero. The values at which this occurs are the initial conditions of the system of figure 22. 
The solution is extracted with the modified Newton method used in the trim procedure. 

With the initial conditions available, the entire system of figure 22 is integrated according to the flow 
diagram in figure 24. The S block is evaluated to get its output y ^ using Ug, which is not necessarily 
the the correct value at the given time. However, the output of S is specified to be independent of 
u s , so the output value is correct. With this input to block C, the values on the right-hand side of the 
figure can be evaluated from the given information on the left-hand side of the figure. With this setup, 
and the input up defined by the user, the integration is carried out using the predictor-corrector method 
in the subroutine DVERK from the IMSL package (ref. 17). A summary of the integration inputs is 
shown in table 6. As in the trim, the variables to be integrated are selected through indices. The input 
is Uj> = Upo + A Up, where Upo is determined in the trim. The term A Up is the input perturbation 
about this trim. 

Linearization 

After trim, it is also possible to linearize a subset of the helicopter model (e.g., block G of fig. 22) 
in the nonrotating system (ISNROT = 1, IINROT = 1, and IONROT = 1). Linearization is done about 
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(225) 


the trim solution using the two-point finite-difference formula, 

df _ L ( x i + Ax j) ~ L ~ Ax J') 

dxj 2 A xj 

The partial-derivative matrix of the helicopter model is defined in figure 25; the rows and columns are 
selected in the same way as in the integration process. For any given time, the partial-derivative matrix 
can be calculated and the linear system can be evaluated according to the following equation: 
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The harmonics of the linear model matrix can be found using 
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These integrals are evaluated numerically using a Gaussian quadrature procedure. The periodic model 
matrix is then expressible as 
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Since nonrotating coordinates are used, the helicopter model matrix will be constant when the model is 
in hover. In forward flight, the model matrix has periodic coefficients, which become more important as 
the forward flight speed is increased. The fundamental frequency u>f depends on the number of blades 
in the helicopter system. If is odd then ujf = n 5 ^, whereas if it is even, uj = A summary of 
the linearization procedure inputs is given in table 7. 


RESULTS 


The model results are compared with other experimental and analytical results as a preliminary 
validation. The first comparison is with blade-equilibrium responses in trim, of a flap-lag-torsion rotor 
model used in reference 8 . Unlike the model described in the present report, the equations in reference 8 
were symbolically generated out to their final form. An omission in the reverse-flow region was 
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discovered in the model of reference 8: the model did not account for the sign change in x ^ in its 
aerodynamic formulation (see fig. 11). This error is introduced into the model described herein so the 
comparison would be valid. The parameters of the model, presented in table 8, are of a four-bladed, 
soft in-plane rotor with a lead-lag frequency of 0.7/rev. The rotor trim is calculated with six harmonics 
( n h T R\i = 6) in each blade degree of freedom, and the fuselage mass is chosen to give a weight 
coefficient of 0.005. The flap, lag, and torsion equations are satisfied for the blade equilibrium while 
simultaneously satisfying the roll-moment, pitch-moment, x-force, and 2 -force equations. The trimmed 
blade-equilibrium response in flap, lag, and torsion is shown in figure 26. Three response curves are 
presented for each blade degree of freedom: (1) a curve from the model from reference 8, which has 
an erroneous reverse-flow formulation; (2) a curve from the model of the present report, which uses the 
same reverse-flow formulation as that used in reference 8; and (3) a curve from the model of the present 
report, which uses what is believed to be the correct interpretation of the aerodynamic formulation 
in the reverse-flow region. The two models with the incorrect reverse-flow interpretation show close 
agreement, which serves to verify the analysis. The alternative reverse-flow interpretation greatly affects 
the torsion response of the system, although its effect on the flap and lag responses is still very small. 

The results from the linearized helicopter model are compared with experimental hover results from 
reference 9. The experimental setup consisted of a three-bladed rotor mounted on a mast that allowed 
pitch and roll motions. Each blade had metal plates with notches, which acted as concentrated hinge 
springs to closely simulate the flap-lag rigid-blade approximation. Torsion springs on the body mount of 
the experimental rotor were adjusted to simulate an air-resonance condition. The parameters necessary 
to model this system are given in table 9. Measurements were made of the damping and frequencies 
of the system modes, and it is to these data that a correlation is made. Figure 27 shows the modal 
frequencies and the lead-lag regressing damping plotted against the rotor rotation rate, with zero blade 
pitch. The simulated air-resonance condition occurs when the body roll-mode frequency nears that of 
the lead-lag regressing mode, which is near a rotor rotational rate of 700 rpm. Excellent correlation with 
the experimental data is seen in this condition. The body pitch and roll damping at zero blade pitch is 
shown in figure 28. Good correlation between experiment and theory is also seen here, although the 
model tends to overpredict the damping at the higher rotor rates. The lead-lag regressing damping at 
9° blade pitch is shown in figure 29; again, the model does a good job of predicting the damping in the 
resonance condition. Detailed cross-plots of the lead-lag regressing damping with respect to the blade 
pitch are shown in figure 30. Good agreement is seen between experiment and theory, although the 
damping tends to fall short at high and negative blade-pitch angles at the higher rotor rates. Figure 31 
shows cross-plots of the body damping at 650 rpm; good agreement in pitch is seen at positive blade- 
pitch angles. The body roll damping tends to be overpredicted, getting worse at higher blade-pitch 
angles. 

The main-rotor wake model used to calculate the downwash on the empennage surfaces is tested 
by correlating it with experimental results from reference 5, page 77. A correlation is done in the 
reference and is also done in the present report by using the flat-wake model described in the Main 
Rotor Downwash subsection. Unlike the model used in reference 5, the model described herein assumes 
a circulation distribution. In reference 5, a blade-element approach is used to find the circulation 
distribution. The integral in equation (152) is evaluated by using a Gaussian quadrature integration 
method along with linearly interpolated values from the pretabulated functions. Figure 32 shows the 
vertical nondimensional downwash at various azimuths at OAR below the rotor. The experimental rotor, 
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with 10° blade twist, is immersed in edgewise flow at Cj — 0.006. Using the flat-wake model gives 
good correlation results; the asymmetry of the downwash at the 90° and 270° azimuth positions is 
captured. 

With the rotor and downwash models validated, it is now possible to validate the full helicopter 
model. The model is configured to simulate a UH-60A helicopter in order to compare model responses 
with the flight-test data that are available in references 10 and 18. Reference 10 presents a validation 
of the real-time model presented in reference 11, which is also known as the GEN HEL model. Ref- 
erence 18 presents a validation of an updated version of the GEN HEL model, which is referred to as 
the Ames GEN HEL model. In the model of the UH-60A used in this report, all automatic controls 
are off (SAS, PBA, etc.). The fuselage aerodynamic representation in the model was generated from 
curve-fitted equations of the static aerodynamic characteristics. The exact equations are given in ref- 
erence 19, pages 2-3, and are used in place of the equivalent-drag-area fuselage aerodynamic model. 
The lag damper is modeled according to reference 11, page 5.1.10, and this representation is substituted 
for the spring/damper model. The control mixer is modeled according to reference 11, pages 5.5.14 to 
5.5.17, without the servo models. The tail-rotor parameters are extracted from reference 20, where a 
detailed description of the UH-60A tail rotor is given. References 11 and 12 are used to extract other 
physical parameters of the helicopter. Table 10 summarizes of all of the inputs to the model. The data 
followed by a question mark are specified for each of the following correlations. From table 10 it can 
be seen that the model has a flap-and-lag main rotor and a six-degree-of-freedom fuselage at constant 
rpm. The dynamic inflow model is included, and the collective pitch, sine pitch, cosine pitch, and 
tail-rotor collective pitch act as the inputs to the helicopter model. No feedback outputs are selected, 
and accelerometer outputs are taken from the fuselage at the positions corresponding to those on the 
flight-test vehicle. 

Trim results calculated with and without the effect of the main-rotor wake in straight and level 
flight are compared with experimental results in figure 33. The trim data were calculated at 25-knot 
increments with the horizontal tail surface at the following settings: 40° for 0, 25, and 50 knots; 11.3° 
for 75 knots; and 4° for 100, 125, and 150 knots. The vehicle was trimmed by satisfying all six fuselage 
force and moment equations, the flap and lag rotor equations, and the inflow equations (see table 10). 
Between 0 and 50 knots, the angle 0 2 was fixed, while the angle 8 X was adjusted (i.e., if = 1,2 in 
table 10). Between 50 and 150 knots, the angle 9 X was fixed, while the angle 6 Z was adjusted (i.e., 
if = 2, 3 in table 10). The collective, longitudinal, lateral, and pedal inputs are presented in percent 
of total motion; the total motion is 10 in. for the stick inputs and 5.38 in. for the pedals. The model 
tends to deviate significantly from the experimental data above 120 knots, where aerodynamic effects 
such as dynamic stall and radial flow, which are not accounted for in this model, can become important. 
Below 120 knots, good correlation is seen and the results are similar to those presented in reference 18, 
p. 28. The lateral and longitudinal stick positions are predicted more accurately by the Ames GEN HEL 
model, but the pedal position is closer for the model described herein. Overall, the results from the Ames 
GEN HEL model are slightly better, probably because that model used tabular data for blade-section 
aerodynamics, as well as other refinements particular to the UH-60A configuration. The flat-wake model 
enhances correlation in all but the collective and lateral stick positions. The collective stick position is 
hardly changed, whereas the lateral stick position is degraded by 1/4 in. at 150 knots. 

Flight test time-histories from reference 10 are compared with the responses of the UH-60A model, 
at 1 knot (figs. 34-37) and at 100 knots (figs. 38-41). Four sets of time-responses are presented at each 
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flight speed; they correspond to 1 in. adjustments to the collective, longitudinal, lateral, and pedal 
inputs. The inputs to the model are those from the flight data, referenced about the trim of the model. 
The trim values generated by the model are indicated by points c, e, a, and p, corresponding to the 
collective, longitudinal, lateral, and pedal inputs, respectively. The angular rates of the fuselage (roll, 
pitch, and yaw) and the accelerometer outputs from the fuselage (longitudinal, lateral, and vertical) 
are also presented for each case. The accelerometer position corresponds to the approximate center of 
gravity of the helicopter; the precise center-of-gravity location is given as POSFUS in table 10. Each 
plot has three curves: the test data, the model response, and the linearized model response. The test 
data are from reference 10 and the test numbers are identical to the numbers used in the reference. The 
model responses come from integration of the full nonlinear system, and the linearized model responses 
come from integration of the linearized model. The responses of the linearized model are added to the 
trim values to give the total responses presented. Since the helicopter is in straight and level flight, 
the angular rates are zero at trim, and the linearized responses are identical to the total responses. The 
accelerometer responses contain gravity effects, and the nonzero trim values of the accelerometers must 
be added to the linearized accelerometer outputs. 

At 1 knot (figs. 34-37), the on-axis angular rates agree well with the test results. The off-axis 
responses of the model do not correlate as well, with the pitch axis being the most troublesome by 
diverging from the test data in tests 212, 203, and 209. However, this is no worse than the results 
from the GEN HEL model, which is refined specifically to the UH-60A configuration and includes an 
engine model with a rotor rpm degree of freedom and actuator dynamics. At 100 knots (figs. 38-41) the 
on-axis responses have fair agreement with experimental results, and again, these results are consistent 
with results from the GEN HEL model. Other tests were run at 60 knots and 140 knots, and the results 
show a strong similarity to the results from the GEN HEL model. Another source for comparison is 
reference 16, in which the results of a linearized model are compared with the same flight data. The 
results from reference 16 are also similar to results obtained with the model of this report. Taking all 
these results together, it is concluded that the present model gives results that are consistent with those 
obtained from the GEN HEL program and the linearized model of reference 16. 

A final comparison is made of the linearized model with the eight-state quasi-static model generated 
from the Ames GEN HEL program. The linearized GEN HEL model is generated with Q = 27 rad/sec, 
ST m = ST C = ST a = 29.583 ft, WL m — WL C = WL a = 20.683 ft, Mp — 492.13 slugs, i a ero = 0, 
and cry = [0,0,0] rad. The linearized model of this report is generated and arranged according to 



(231) 

(232) 


The vector x r contains all rotor states and inflow states, and the vector xj contains all fuselage states. 
The inputs are the longitudinal, collective, lateral, and pedal positions (uj = [S e , 6 C , 8 a ,8 p ) T ). The out- 
puts y are chosen to correspond to the states of the GEN HEL model (y r = [u, w , q, Ad, v, p, A0, r, At/?]). 
For comparison with the GEN HEL model, the derivatives of the rotor and inflow states x r are set to 
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zero, yielding 


x_f Ay^^^Hf 

(233) 

V ~ C re dS.f D re dU 

(234) 

^red = Aff ” Aj r A rr A r f 

(235) 

Bred ~ B j — Af r A rr Br 

(236) 

C 'red = C f — C t A tt A r f 

(237) 

Bred ~ B — CrA rr Br 

(238) 


The state variables used to describe the fuselage motions are not precisely those used in the GEN HEL 
model; however, the information needed to make the transformation to the GEN HEL states is found in 
the equation 
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' d x ' 
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0 Z 
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> — [C* re d] < 
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Ox 
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k 0z , 


(239) 


The quantities u, v, and w are the velocities along the fuselage, and p, q, and r are the body angular 
rates. The quantities A0, A <f>, and Aip are the pitch, roll, and yaw angular perturbations. With this 
choice of output quantities, it is clear that D re d = 0, and thus equation (239) holds. Equation (239) is 
then used to transform the model to the final form 


' u 


u 

w 


w 

Q 


Q 

AO 

> = F< 

a e 

V 


V 

V 


V 

A<\> 


A (f> 

, r j 


, r , 


+ GI 


be 

b c 

ba 

h 


> 


(240) 


where A tp has been removed from the system. 


This operation was carried out on a linearized model at l knot; the resulting F and G matrices 
are presented in table 11 along with the system matrices from the Ames GEN HEL program. The 
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eigenvalues and normalized eigenvectors are shown in figure 42 for both models. The calculations are 
based on the F matrix represented in dimensions of feet, seconds and degrees so as to scale the angular 
rates to an order of magnitude comparable to the velocities. Real eigenvalues have real eigenvectors 
(figs. 42(d) and 42(e)); complex eigenvalues have complex eigenvectors composed of a real part and 
an imaginary part. The eigenvector corresponding to the complex conjugate of the eigenvalue is the 
complex conjugate of the eigenvector. The greatest difference between models arises in the higher 
frequency modes, shown in figures 42(a) and 42(b). The model of this report tends to introduce more 
coupling between roll and pitch than the Ames GEN HEL model. The mode in figure 42(c) is dominated 
by w and r in both models, although the content differs. In figures 42(d) and 42(e) the two modes are 
very close, although the damping is much higher in the model of this report. Differences in the damping 
and frequencies of the eigenvalues are also noted; they are likely caused by the different inflow models 
used in the two analyses. Selected Bode plots of the system are shown in figure 43. The 8 c /w plot 
shows good agreement and the 6 p /r gain agrees well with the Ames GEN HEL model. The 360° phase 
shift in the 8 p /r plot is due to a zero in the right half of the s-plane that is not unstable in the model of 
this report. The sharp peak in the 6 a /p plot is due to the light damping that is predicted by the Ames 
GEN HEL model. 

The entire process is repeated at 100 knots, resulting in the F and G matrices given in table 12. The 
eigenvalues and eigenvectors are presented in figure 44; the eigenvectors agree more closely between 
models in this case than in the 1-knot case. The on-axis Bode plots of the 100-knot case are given in 
figure 45. The two models do not agree closely, but the general trends in the models are similar. 

CONCLUDING REMARKS 

A single-rotor helicopter mathematical model for the design of high-bandwidth flight control sys- 
tems has been described. The flap, lag, and torsion rotor model allows the approximate modeling of 
hingeless rotor systems, and the various hinge sequences that are available allow the modeling of many 
different articulated rotors. Quasi-steady Greenberg aerodynamics are used in conjunction with the 
nonlinear Pitt/Peters dynamic inflow model to calculate the blade aerodynamic forces. An rpm degree 
of freedom of the main rotor is available, along with a six-degree-of-freedom rigid fuselage. The tail 
surfaces have main-rotor downwash effects that arise from the flat-wake simulation. Operations on the 
model include trim, linearization, and time-integration, which have been formulated so that they may 
be applied to any consistently cast rotorcraft mathematical model. This features allows for modification 
of the helicopter model without requiring significant changes in the operations. Preliminary validation 
of the model showed reasonable correlation with experimental and analytical results. 

Further validation of the model is necessary. The higher order linearized models need to be 
compared to experimental results, and flap-lag-torsion experimental data are needed to validate the 
torsional effects of the rotor model. Also, additional work is necessary to extend the range of valid flight 
conditions, which are currently limited to low-speed, lightly loaded helicopter systems. The limiting 
factor is the blade-section aerodynamic model, which should be improved to account for dynamic 
stall, radial-flow effects, and compressibility effects. Finally, the linearization process is limited to the 
helicopter model itself. A useful additional operation would be the linearization of the remaining blocks 
of the system and the combination of these linearized blocks into a single set of system matrices. 
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Table 1. Rotation sequences of blade coordinate systems 
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P to L 

L to B 
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Table 2. Function odd/even properties 


Function 

Even/odd (-r i , yi - ~ i) 

Equation number 

Figure number 

-\i'x 

(E,0,E) 

2.10 

2.15 

Xvy 

(0,E,E) 

2.9 

2.14 

M 

(0,E,0) 

2.37 

2.25 

X 

(0,E,E) 

2.40 

2.26 

L 

(0,0,0) 

2.38 

2.28 

K 

(0,0,E) 

2.39 

2.27 
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Table T Summary of helicopter propertic 



38 




39 



Tabic.’}. Continued 



Read in flat wake table for Tit 
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Symbol Computer Description Dimensions 



nondimensional, M = mass, L = length, rad = radians, 1' = time. 



Table 1- Information passed to driver from model 



X 

5b 
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Summary of trim variables- namelist TRMDAT 



I 


.linearization 



Table (). Summary of integration variables namelist INTDA'l 
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Table 7. Summary of linearization variables namelist LIN DAT 



coefficient matrices 

Nil LIN Number of harmonics in periodic system linear coefficient matrices 

IOLYVRT Write out coefficient matrices in readable format? (0 — no, I = ye? 
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Value 

Rotor/fuselage 
3 
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Table 10. Continued 
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Table 10. Concluded 
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Table 1 l. System matrices from CIEN 1 1 EL program and model of this report at l knot 
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Table 12. System matrices from GEN HEL program and model of this report at 100 knots 
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Figure 1. Fuselage-fixed coordinate frame (F). 
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Figure 2. Free-body diagram of fuselage. 
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Figure 4. Fuselage wind-axis coordinate frame (W) and its transformation to the F frame. 
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Figure 7. Flap-lag-pitch-torsion hinge sequence. 
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Figure 8. Lag-flap-pitch-torsion hinge sequence. 
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Figure II. Blade-section aerodynamics. 













Figure 14. Tail-surface drag coefficient break points. 
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Figure 16. Tail-rotor orientation and aerodynamic forces. 
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Figure 21. General input/output model structure used by trim, linearization, and integration operations. 
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Figure 22. General system structure for integration 
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Figure 25. Partial-derivative matrix for model linearization 
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Figure 26. Blade responses from the present model compared with those from a flap-lag-torsion rotor 
in trim. 
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Figure 27. Modal frequencies and lead-lag regressing damping data from the present model compared 
with those from a flap-lag rotor/fuselage system at zero blade pitch. 
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Figure 28. Body pitch and roll damping data from the present model compared with those from a 
flap-lag rotor/fuselage system at zero blade pitch. 
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ROTOR RATE, rpm 

Figure 29. Lead-lag regressing damping data from the present model compared with those from a 
flap-lag rotor/fuselage system at 9° blade pitch. 
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Figure 30. Lead-lag regressing damping data from the present model compared with those from a 
flap-lag rotor/fuselage system at various rotor speeds and blade pitches. 
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Figure 31. Body pitch and roll damping data from the present model compared with those from 
flap-lag rotor/fuselage system at 650 rpm. 
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FLIGHT SPEED, knots FLIGHT SPEED, knots 

Figure 33. Trim in straight and level flight of helicopter model compared with llight-test results. 
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Figure 36. UH-60A model responses with l-in.-left lateral stick at 1 knot. 
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Figure 37. Concluded. 
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Figure 38. Concluded. 






YAW RATE, deg/sec PITCH RATE, deg/sec ROLL RATE, deg/sac 































VERT ACC. ft/sec* LAT ACC, ft/sec 2 LONG ACC, ft/sec 


AEFA TEST 102 


(- 


u 


0 



LINEARIZED f ~ . ^ 

MODEL 




0 1 2 3 4 5 6 


TIME, sec 


il = 27.023 rad/sec M p = 465.81 slugs ST m = ST C = ST a = 29.058 ft f(4) 


'f= 


Aero 


4659 

0 

0 

= 0 


0 

38512 

-1882 


0 

-1882 

38512 


WL m =WL c =WL a =19.29 ft i f 
Oy= 0.0698, 0.0698, 0 rad 


IOPT5 = 5,5,5 


LINEARIZED, 

MODEL 


161.6 ft/sec 
2,3 


Figure 41. Concluded. 
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Figure 42. Eigenvalues and eigenvectors of this report’s model and the Ames GEN HEL model at 
1 knot. 
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Figure 44. Eigenvalues and eigenvectors of present model and Ames GEN HEL model at 100 knots. 
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